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ABSTRACT 


Maxvvorthy's  experiments  with  the  flight  of  a  sphere  through  an  incom¬ 
pressible  electrically  conducting  fluid  with  an  aligned  magnetic  field  showed 
that  transition  in  the  trailing  wake  was  moved  to  higher  fluid  mechanical 
Reynolds  numbers.  Theoretical  and  experimental  work  on  suppression  of 
vorticity  in  uniformly  conducting  incompressible  fluids  with  imposed  uniform 
transverse  magnetic  fields  is  analyzed  to  find  a  correlating  parameter  for  Cq, 

Cf  and  transition  to  turbulence.  The  parameter  which  is  found  is  the  square 
root  of  the  ratio  of  the  fluid  mechanical  Reynolds  number  (R)  to  the  electromag¬ 
netic  interaction  parameter  (S),  n/r/S,  for  the  conditions  S«  0(1)  or  RS»0(1). 

The  is  the  Reynolds  number  based  on  the  Hartmann  layer  thickness, 

R.I.J  =  sj p V 2 /ctB 2 v\  It  is  shown  that  for  the  cases  RS>>0(1),  the  absolute  value  of 
Rpj  for  transition  is  of  the  same  order  as  the  absolute  value  of  R  for  transition 
in  the  zerofield  case,  Rx  0  •  Thus,  the  apparent  effect  of  the  transverse  mag¬ 
netic  field  in  the  incompressible  case  is  to  move  Cj),  Cf  and  transition  to  higher 
values  of  the  fluid  mechanical  Reynolds  number,  Rg,  when  the  field  is  on.  The 
apparent  correlation  leads  to  Rg  =  RqIRq  ‘ 

Extension  of  the  incompressible  problem  to  the  compressible  one  is 
based  on  the  observed  similarity  of  the  transition  processes  in  the  hypersonic 
and  incompressible  wakes  as  suggested  by  Fay  and  Goldburg.  Average  values 
of  the  fluid  properties  in  the  base  region  were  chosen  as  typical  for  the  calcula¬ 
tion  of  hypersonic  wake  transition  Reynolds  numbers,  R  and  Rpj.  The  presence 
of  large  amounts  of  easily  ionizable  ablation  products  were  assumed  to  give 
Spitzer  conductivity.  The  calculated  values  for  the  relevant  nondimensional  fluid 
mechanical  and  electromagnetic  parameters  for  typical  ballistic  range  conditions 
indicated  they  were  of  requisite  magnitude  to  affect  the  flow. 

Hypersonic  experiments  were  performed  covering  the  range  from  R  ^5Rx  q 
with  zero  field  to  Rh  ^/^R-X  0’  where  Rx,  0  is  the  zero  field  wake  transition 
Reynolds  number.  The  magneto-hypersonic  flow  fields  were  produced  by  0.  22  cali¬ 
ber  nylon  pellets  traveling  at  about  14,000  feet  per  second  iu  argon  through  the  axial 
magnetic  field  of  a  solenoid  produced  by  the  discharge  of  a  capacitor  bank.  Drum 
camera  pictures  of  the  self-luminous  wakes  and  photomultiplier  traces  were  used 
as  the  primary  diagnostic  tools--the  self-luminous  ablation  products  color  the  wake 
in  the  same  fashion  as  dye.  The  camera  was  oriented  such  that  its  field  of  view  en¬ 
compassed  the  region  before  and  including  part  of  the  coil. 

In  the  first  set  of  runs  made  at  5  cm  Hg  pressure,  a  change  in  the  lumin¬ 
ous  structure  of  the  flow  field  occurred  possibly  suggestive  of  a  change  from 
turbulent  to  laminar  wake  flow. 

In  the  second  set  of  experiments  made  at  5  cm  and  2  cm  Hg  pressure,  no 
MHD  effect  on  the  hypersonic  turbulent  flow  field  was  evident.  The  paper  con¬ 
cludes  with  a  discussion  of  these  results. 


-Ill- 


I.  Correlation  Parameter  for  MHD  Flows  with  Transverse  Fields 
A.  Introduction 

The  macroscopic  magnetohydrodynamic  equation  of  motion  for  an 
electrically  conducting  incompressible  medium  in  which  the  Hall  current 
effect  and  the  ion  slip  effect  are  neglected  and  in  which  the  current  loops 
close  on  themselves  is 

DV  1  1  2 

-  +  -  grad  p  = —  _j  X  ^  V  V 

dt  p  p 

2  =  a(E  +  V  X  B) 

in  the  usual  notation.^  For  scalar  conductivity  the_2  x  B  force  is  seen  to 
have  a  component  directed  oppositely  to  V,  of  order  cTVB^. 

.  In  shear  flow,  ay  shown  in  Fig.  1,  a  uniform  transverse  field  with 
a  uniform  conducting  medium  will  cause  eddy  currents  to  circulate  such 
that  the  fast  fluid  decelerates  more  strongly  than  the  slow.  The  result  is 
to  suppress  the  vorticity  in  the  flow  field.  In  this  argument,  the  assump¬ 
tion  is  made  that  the  induced  field,  B^,  caused  by  the  eddy  currents  is 
small  relative  to  the  imposed  field  Bq-  Figure  2  shows  a  schematic  of  a 
B  field  interaction  with  the  vortex  ring  in  the  base  of  a  body.  It  is  seen 
that  the  component  of  the_2  x  B  force  oppositely  directed  to  V  slows  down 
the  speed  of  rotation  of  the  ring  cross  section  causing  a  decrease  in  the 
strength  of  the  base  vortex  ring. 

The  MHD  equation  of  motion  (1)  is  nondimensionalized  such  that 
each  term  of  the  equation  is  composed  of  two  factors:  a  nondimensional 
parameter  indicating  the  term's  physical  significance  and  its  magnitude 
relative  to  the  other  terms,  and  a  function  of  the  unknowns  and  their  de¬ 
rivatives  of  magnitude  of  order  one.  This  nondimensionalization  results 
in  the  following: 

-;3-  +  grad  p  =  —  .  V  V  +  S  •  j  x  B  (2) 

dt  R 


where  the  tilde  represents  the  usual  nondimensional  quantities  and 


R  = 


V  C 
o  o 


..u  r  f  {  inertia  | 

represents  the  ratio  of  forces  < — s -  / 

^  I  VISCOUS) 

per  unit  volume 


(3) 


S  = 


%®o  ^o  ) 

- ^ -  represents  the  ratio  of  forces  <• 

P  r-\  ' 

per  unit  volume 


o  o 


j  X  B 
inertia 


(4) 


The  single  subscript  (o)  represents  a  characteristic  quantity.  The  charac¬ 
teristic  length  governing  a  problem,  Hq,  is  the  length  over  which  the  forces 
interact  to  cause  the  significant  change  in  velocity.  In  deducing  S,  the  as¬ 
sumption  has  been  made  that  the  velocity  perpendicular  to  Bq  is  like  V^- 

If  1/R  and  S  are  large,  then  the  inertia,  terms  on  the  left-hand  side 
of  Eq.  (Z)  can  be  neglected  and  another  nondimensional  parameter,  (R-  S), 
appears;  the  \!  R  •  S  is  usually  called  the  Hartmann  number, 


o  B 
o  o 


V 

o 


represents  the  ratio  of  forces 
per  unit  volume 


j  ^  \ 

viscous  ) 


(5) 


The  magnetic  Reynolds  number,  Rj^l-  is  the  nondimensional  paramete  r 
which  represents  the  severity  of  the  effect  of  the  electrically  conducting 
fluid  flow  on  the  magnetic  field. 


=  ^o%^o^o  represents  j^|  (6) 


where  Pq  is  the  fundamental  constant,  permeability  of  vacuum,  in  MKS  units. 

In  a  real  gas  the  electrons  tend  to  spiral  aroui  d  magnetic  field  lines 
between  collisions  giving  rise  to  a  tensor  conductivity  rotating  and  decreas¬ 
ing  the  current  vector  and  hence  the  x  B  force  vector.  The  Hall  coefficient, 
tor  >  is  a  measure  of  the  angle  between  the  V  x  B  vector  and  the  current  vec¬ 
tor  due  to  the  tensor  conductivity. 


U)T 


o  B 


o  o 


N  •  e 
e 


represents  ratio  of  times 


period  for  collision  ) 
period  for  revolution / 
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where  ou  is  the  electron  spiraling  frequency  along  field  lines  and  j  is  the 
mean  free  time  between  collisions.  Ng  is  the  electron  concentration  and 
e  the  electron  charge.  In  liquid  metals,  this  term  is  effectively  zero. 

To  effect  an  MHD  vorticity  interaction  the  magnitudes  of  the  non- 
dimensional  parameters  should  be  of  the  following  order 


0(1); 


<  0(1);  and  OJT  <  0(1). 


It  may  be  that  the  condition,  the  ratio  of  |  j  x  B/inertia|  forces  be  of  order 
one  [S  ^  0(1)],  can  be  replaced  by  the  less  severe  condition,  the  ratio  of 
I  j  X  B/viscous  I  forces  be  greater  than  of  order  one;  RS  =  ^  0(  !)• 

Available  experiments^'^'  ^  do  not  encompass  .a  sufficient  range  of  param¬ 
eter  magnitudes  to  answer  this  question  definitely.  However,  these  exper¬ 
iments  do  show  some  evidence  that  there  is  an  effect  when  S  is  less  than 
order  one,  so  long  as  RS  »  0(1). 

B .  Theory 

Under  the  conditions  of  steady  flow,  with  a  homogeneous  B  field  in 
only  the  y  direction,  Lundquist"^  following  a  scheme  of  Lehnert^  was  able 
to  reduce  by  linearization  the  comp’  te  set  of  equations  governing  MHD 
phenomena  to  the  following  "boundary  layer"  form; 


2 

8 


U  O’ 
^o  o 


V 


The  first  term  represents  the  electromagnetic  driving  force  and  the  second 
term  represents  the  electromechanical  inte ’•action.  This  "boundary  region" 
equation  indicates  that  changes  in  V  along  tiie  magnetic  field  direction  take 
place  in  length  of  order  Lj^l. 


Lundquist  then  surmises;  The  fact  that  the  variation  of  the  velocity 
is  confined  to  the  "boundary  region"  thickness,  suggests  that  the  length 

dimension  entering  the  ordinary  fluid  mechanical  Reynolds  number  be  re - 
placed  by  Thus,  the  correlating  parameter  in  the  MHD  case,  R^, 


becomes: 


V 


R 


H 


L 


V  (e  ) 
o  o 


P  V 

o 


M 


o  B  ) 

o  o  o 


=  ^yR/s 


(9) 


where  now  R  and  S  are  based  on  any  characteristic  length,  say  a  body  diam¬ 
eter  or  channel  diameter,  D.  The  nondimensi  jnal  parameter,  \l  R/S , 
represented  by  the  symbol,  Rh>  is  the  fluid  mechanical  Reynolds  number 
based  on  the  Hartmann  layer  thickness  Ljvl- 

One  can  look  at  this  another  way.  In  the  MHD  equation  of  motion  (2) 
there  are  two  independent  nondimensional  parameters,  R  and  S.  These 
can  be  rearranged  into  =  RS  and  either  R  or  S.  It  can  be  claimed 

that  the  interactions  of  interest  take  place  over  a  region  in  which 

=  I  j  X  B/viscous  [  =  0(1).  This  region  is  of  extent  Lj^,  thus  recover¬ 
ing  Eqs.  (8)  and  (9). 

'  Lykoudis^  noticed  the  similarity  of  the  velocity  profiles  in  the  case 

of  the  laminar  sublayer  with  mass  suction  and  the  laminar  channel  flow 
v’ith  transverse  magnetic  field.  By  taking  over  the  results  for  the  termin¬ 
ation  of  the  laminar  sublayer,  he  was  able  to  establish  a  criterion  for  trans¬ 
ition  in  the  MHD  channel  flow  case.  This  criterion  depended  solely  on  the 
parameter  R/HjvI  ='s/  R/S  =  RH-  And,  by  using  the  value  of  the  transition 
Reynolds  number  for  the  zero  field  laminar  sublayer,  he  was  able  to  pre¬ 
dict  the  value  of  RHjX>  which  agreed  with  Murgatroyd's  experiments  to  be 
described  next.  ' 

C.  Incompressible  Experiment 

The  available  relevant  experiments  are  incompressible:  (1)  channel 
flow  with  transverse  magnetic  field  and  (2)  a  sphere  in  free  flight  with  mag¬ 
netic  field  aligned  with  the  flight  direction.  In  this  discussion  S  is  based  on 
the  assumption  that  the  velocity  perpendicular  to  B  is  of  the  order  of  the 
free  stream  velocity.  This  is  satisfied  in  the  channel  flow  with  transverse 
field,  and  also  in  the  sphere  base  flow  with  field  aligned  to  the  flight  direc¬ 
tion  since  the  base  vortex  ring  has  velocities  trans/erse  to  the  field  which 
are  of  the  order  of  the  free  stream  velocities.  Fig.  2.  Further  downstream 
behind  the  sphere  the  flow  tends  to  be  parallel  to  the  field  direction.  MHD 
vorticity  interaction  with  instabilities  in  the  laminar  far  wake  with  a  flight 
aligned  magnetic  field  is  a  separate  problem  which  would  have  as  its  chan¬ 
nel  analogue  pipe  flow  transition  with  a  flow  aligned  magnetic  field.  (This 
last  problem  has  been  attacked  experimentally  by  Globe  '  and  theoretically 
by  Stuart;^  and  the  MHD  effect  was  significantly  less  than  that  to  be  de¬ 
rived  below.  ) 
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Murgatroyd^  performed  experiments  in  which  measurements  of 
pressure  gradients  were  taken  during  the  flow  of  mercury  in  a  rectangular 
channel  of  15  :  1  aspect  ratio  in  the  presence  of  a  transverse  magnetic 
field.  The  experiment  was  designed  to  study  the  relationship  between 
pressure  gradient,  mean  flow  velocity,  and  field  strength.  He  found  that 
it  was  possible  to  suppress  turbult  ■•'.ce  even  at  the  highest  obtainable  fluid 
mechanical  Reynolds  number  (R  ~  25,000  where  the  characteristic  length 
was  taken  as  the  channel  half -width). 

Murgatroyd  referenced  Lundquist' s  reasoning  given  in  the  previous 
section  and  plotted  his  data  as  Cf,  the  usual  skin  friction  coefficient,  agai 
Hj^/R  =  (^y  R/S)“  i,  Fig.  3.  These  experiments  show  that  over  a  range  of 
from  1  to  0.  4  and  for_R  based  on  channel  half  width  of  order  25,  000,  is 
a  function  of  Rjq  -  R/S  alone  for  \/ R  •  S  >  100.  Transition  was  noted  as 
occurring  at  a  definite  value  of  the  parameter,  Rh)x  ~  N/'  R/S  ^  ~  225. 

like  25,000.  This  agrees  with  Lundquist' s  criterion,  since  for 
the  ordinary  fluid  mechanical  zero  field  case,  R)x,  B=  0  200.  (The 

value  of  Rh))^  predicted  by  Lydoudis^  was  236.  ) 

Maxworthy^  performed  a  series  of  experiments  in  which  the  drag 
experienced  by  metal  spheres  of  several  different  diameters  was  measured 
by  determining  their  terminal  velocities  as  they  fall,  vertically,  through 
an  electrically  conducting  fluid  (liquid  sodium)  and  a  solenoid  magnetic 
•field.  Induction  coils  were  used  to  detect  the  moving,  perturbed  magne'ic 
field  associated  with  the  sphere,  allowing  accurate  determination  of  the 
latter's  position  in  space  and  time  and  the  gross  nature  of  the  perturbed 
field. 


Now,  according  to  Lundquist' s  suggestion  we  expect  the  onset  of 
wake  uns teadine s s  to  occur  in  the  MHD  sphere  wake  at  the  value  of 
RH  =  ^  R/S  equal  to  that  value  of  R  at  which  it  occurs  in  the  zero  field 
case.  The  experiments  of  Magarvey  and  Bishop  as  analyzed  by  Fay  and 
Goldburg9  yield  (R  based  on  diame  te  r) 


R 


crit) 


200. 


B  =  0 


(10) 


Thus ,  the  prediction  in  the  MHD  case  for  S  =  0(1)  or  RS  »  0(  1)  is 


crit 


200. 


(11) 


Maxworthy' s  experiments  tend  to  confirm  this.  Table  II  is  a  list¬ 
ing  of  the  test  conditions  for  11  flights.  The  appears  to  be  between 

111  and  143.  At  Rh  ~  100  the  coil  output  is  smooth  and  at  Rh  >  200  the 
coil  output  shows  pe  riodicity.  2  The  standard  quantitative  measurement 


-5- 


c  in 


of  wake  periodicity  in  the  zero  field  case  is  the  behavior  of  the  Strouhal 
number  (Sj)  defined  at  the  "shedding"  frequency  (f)  times  the  body  diam¬ 
eter  (D)  divided  by  the  free  stream  velocity  (V)  with  R.  The  behavior 
for  the  zero  field  case  is  discussed  in  Fay  and  Goldburg9  where  it  is  shown 
that  a  Strouhal  number  of  0.  1  is  observed  at  an  R  of  about  300.  Ma.xworthy 
points  out  that  the  Strouhal  number  of  the  run  at  an  Rh  210  is  also  0.  1. 

It  is  important  to  note  that  the  fluid  mechanical  diameter  Reynolds  number, 
R,in  all  of  these  cases  is  greater  than  an  order  of  magnitude  above  Rcrit 
for  the  B  =  0  case:  5  to  10  thousand  compared  with  200. 


Table  II.  MHD  SPHERE  WAKE^ 


,  =  \/R/S 

Dia. 

Downstream 

W  ake 

S 

27.  8 

5330 

Steady 

6.  95 

35.  7 

2140 

Steady 

1.  68 

62.  5 

7270 

Steady 

1.  60 

105.  0 

‘l  1 1 .  0 

5600 

Steady 

Steady 

0.  48 

143.  0 

3290 

? 

0.  161 

143.  0 

10000 

Unsteady 

0.  49 

210.  0 

11200 

(Strouhal=0.  1  0) 

0.  24 

222.  0 

11000 

Unsteady 

0.  223 

286.  0 

21000 

Unsteady 

0.  254 

500.  0 

3720 

Unsteady 

0.  015 

R 


crit 


B 


a  200  (Ref.  9). 

0 


Figure  4  shows  Cj)  for  a  sphere  as  a  function  of  R  alone  when  B  =  0 
(the  classical  result^^O),  and  as  a  function  of  Rpj.  For  the  conditions  of  Max- 
worthy's  experiments,  RS;^400,  it  is  seen  that  Rh  correlates  the  data  well 
(parenthetically,  S  <  0(1)  ),  and  that  for  the  same  values  of  Rh  and  R)2  =  o 
the  corresponding  values  of  Cd  are  always  within  a  factor  of  two  of  each 
othe  r. 
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D.  The  Number  Rh  =  n/ R/S 

The  above  evidence  suggests  th^t  an  indication  of  how  the  MHD 
vorticity  interaction  problem  scales  with  an  S  ~  0(1)  or  RS  »  0(1)  may 
be  obtained  by  replacing  the  usual  fluid  mechanical  Reynolds  numbe  r  of 
the  zero  field  case  with  the  Hartmann  layer  Reynolds  number,  Rj^=\/  R/S. 
This  broad  role  for  the  number  Rh  prompts  one  to  search  for  a  general, 
though  perhaps  imprecise,  deduction  of  Rh  from  the  equations  of  motion. 
Equation  (2)  might  be  characterized  in  the  following  way: 


(inertia  terms) 


=  \i 

+  1  s 


viscous  forces 
inertia  forces 

/  j  X  B  forces 
I  inertia  forces 


(viscous  terms) 


(j  X  B  terms) 


In  ordinary  incompressible  viscous  fluid  mechanics,  most  effects  scale 
with  the  ratio  of  the  first  power  of  driving  force  to  retarding  force,  i.  e.  ,  R. 
In  the  combined  MHD  vorticity  interaction  case,  one  can  make  a  simi¬ 
lar  number  in  two  ways:  an  arithmetic  combination,  or  a  geometric  com- 
laination,  where  q  is  an  integer  between  0  and  1.  The  following 

table  lists  three  geometric  combinations. 

Correlating  Parameter 


ordinary  viscous  0  1  R 

combined  viscous  0(1)  1/2 

magneto 

ordinary  inviscid  <x>  0  1 /S 

magneto 


It  is  curious  that  although  this  correlating  parameter  has  appeared 
in  an  important  way  in  both  theoretical  and  experimental  investigations,  it 
has  not  been  singled  out  as  such.  The  q  =  l/2  case  has  u^'ually  been  written 
in  terms  of  the  Hartmann  number  Rh  =  Hm/S.  This  seems  to  obscure  its 
fundamental  nature  as  the  ratio: 


^  _ /  (  reta 


driving  force 
etarding  fore 


Type  of 

Fluid  Dynamics 


/  j  X  B  forces  i 
(viscous  forces! 


inertia 

viscous 


..  “11/2  r 

inertia  ' 

j  X  B  J 


driving  force 
retarding  force 


(13) 
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II. 


Hypersonic  Magnetogas  Dynamic  Wake  Experiments 

In  the  above  discussion  we  have  suggested  that  an  indication  of  how 
the  MHD  vorticity  inte  raction  problem  scales  with  an  S  0(1)  or  RS  >0(1) 
may  be  obtained  by  replacing  the  usual  fluid  mechanical  Reynolds  numbe r 
of  the  zero  field  case  with  the  nondimensional  parameter,  Rpj  =  R/S.  For 
compressible  flow,  it  is  usual  to  extend  the  incompressible  result  to  the 
compressible  case  with  the  qualification  that  the  representative  local  values 
rather  than  free  stream  values  are  used.  Values  characteristic  of  the  base 
region  conditions  are  taken  as  representative  for  the  present  purposes. 

A  hypersonic  counterpart  of  Maxworthy' s  incompressible  experiment 
is  a, ballistic  range  experiment  in  which  pellets  fly  through  an  aligned  sole¬ 
noid  magnetic  field.  Evaluation  of  the  relevant  MHD  parameters  in  the  base 
region  of  hypersonic  pellets  was  made  for  the  following  assumptions  of  rep¬ 
resentative  aerodynamics  base  region  conditions  ('b). 


5P 

00 

(Al) 

V 

00 

(AZ) 

l/Z  h  , 

stagnation 

(A3) 

Equilibrium  Chemistry 

(A4) 

N  = 
e 

l/lO  giving  Spitzer 

conductivity 

(A5) 

where  h  is  the  enthalpy  and  is  the  total  number  of  particles  per  unit  vol¬ 
ume.  Argon  gas  is  taken  to  be  the  working  fluid  to  maximize  gas  luminosity 
and  electrical  conductivity. 

Calculations  for  the  base  region  values  for  the  parameters  S/B^  i, 
Hm/B£  and  Rh  '  B  are  shown  on  a  pressure-velocity  grid  in  Fig.  5; 

{S/£)a;T  =  1.  0.  =  1.0  and  (RH)q;T  =  1.  0  In  Fig.  6,  and  a)T  /B  and 

rm/j2  in  Fig.  7.  (Al)  and  (AZ)  are  justified  on  the  basis  of  theory  and  ex¬ 
periment,  the  velocity  referring  to  the  expansion  of  flow  around  the  flare 
and  along  the  base  region.  (A3)  has  not  been  checked  definitively,  but  it 
seems  reasonable  and  is  a  commonly  used  representative  value.  (A4) 
is  based  on  a  long  residence  time  for  fluid  in  the  base  region.  (A5)  is  the 
most  questionable  assumption.  It  is  based  on  the  assumed  presence  of 
large  amounts  of  easily  ionizable  ablation  products. 

Fay  and  Goldburg9  indicate  that  transition  from  laminar  to  turbulent 
wakes  behind  0.  ZZ-inch  diameter  spheres  in  argon  at  velocities  like  14,000 
ft/sec  occurs  at  a  Reynolds  number  based  on  shoulder  conditions  of  ZOOO. 

The  corresponding  base  region  value  is  like  600.  The  corresponding  ambient 
pressure  is  like  1  cm  Hg. 
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The  suggested  criterion  for  suppression  of  wake  transition  is  that 
the  fluid  mechanical  Reynolds  number  based  on  Hartmann  layer  thickness, 
RH.  be  less  than  the  ordinary  fluid  mechanical  Reynolds  number  R^,  0- 
For  a  range  pressure  of  5  cm  Hg  and  velocitie s  like  14,000  ft/sec,  this 
condition  is  met  when  the  magnetic  field  strength  B  is  chosen  to  make 
OJT  =  1.0:  Rh  «  Rx,o  ~  ^  0(1).  RS  «  0(10^))  fulfilling  the 

requirements  for  an  MHD  vorticity  interaction. 

A  set  of  MHD  parameter  maps  assuming  that  the  gas  in  the  base 
region  is  ionized  argon  alone  is  contained  in  Appendix  A.  These  maps 
indicate  that  an  experiment  in  the  ballistic  range  is  not  feasible  if  the  gas 
properties  are  assumed  to  be  due  solely  to  ionized  argon.  Fig.  A-Z.  In  the 
pure  argon  case,  for  the  same  ballistic  range  conditions,  RpjT  ^>1000. 

The  experiment  to  explore  the  possibility  of  suppressing  hypersonic 
wake  transition  is  shown  in  Fig.  8.  Nylon  pellets  of  0.  22-inch  diameter 
are  fired  at  velocities  near  14,  000  ft/sec  through  the  magnetic  field  of  a 
solenoid  placed  about  a  glass  test  section  in  the  AERL  ballistic  range.  The 
magnetic  field  of  the  solenoid  is  produced  by  passing  a  current  through  the 
wire  by  the  discharge  of  a  capacitor  bank.  The  inductance  of  the  coil  and 
capacitance  of  the  bank  were  chosen  so  that  the  discharge  quarter  cycle 
time  would  be  130  microseconds.  This  ringing  time  of  the  capacitor  bank 
is  sufficiently  long  to  produce  nearly  DC  current,  while  the  pellet  is  in  the 
‘coil.  A  maximum  magnetic  field  strength  of  2.  6  webe  r s/mete  r^  is  obtained 
on  the  coil  centerline  with  the  present  design.  The  test  and  coil  section  are 
made  of  transparent  material,  pyrex  and  Incite,  in  order  that  streak  photo¬ 
graphs  can  be  taken  when  the  pellet  is  inside  of  the  solenoid.  Fig.  9.  The 
use  of  streak  drum  camera  photographs  to  differentiate  between  laminar 
and  turbulent  wakes^’  has  been  reported  previously.  The  streak  camera 
is  oriented  such  that  its  field  of  view  encompasses  a  region  before  and  in¬ 
cluding  part  of  the  coil.  The  appearance  of  luminous  structure  in  drum 
camera  photographs  is  interpreted  as  turbulent  wake  flow  whereas  smooth 
and  diffuse  luminosity  is  interpreted  as  laminar  flow. 

The  base  MHD  experiment  was  composed  of  two  sets  of  runs  sep¬ 
arated  by  runs  for  another  experiment  on  the  AERL  range.  In  the  first 
set  of  runs,  a  single  streak  drum  camera  was  oriented  such  that  its  field 
of  view  encompassed  a  region  before  and  including  approximately  one -half 
of  the  coil.  In  the  second  set,  a  second  drum  camera  was  added  so  that 
all  of  the  coil  region  as  well  as  some  region  upstream  of  the  coil  could  be 
photographed. 

The  very  first  runs  were  made  with  the  coil  unenergized.  A  laminar 
and  turbulent  streak  photograph  is  shown  in  Fig.  10  and  11,  respectively. 
The  streak  photographs  of  the  pellet  within  the  coil  provides  sufficient  in¬ 
formation  to  reveal  the  absence  or  presence  of  luminous  turbulent  structure 
in  the  wake.  The  initial  vertical  break  in  the  streak  photograph  is  the  down¬ 
stream  coil  support  face  plate.  The  subsequent  ve  rtical  b  reaks  are  the  wire 
turns  of  the  solenoid.  The  luminous  streaks  which  characterize  turbulent 
flow  can  be  recognized  within  the  coil  in  Fig.  11. 
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A  streak  photograph  with  the  coil  energized  is  shown  in  Fig.  12. 

Such  a  photograph  was  obtained  in  each  of  eight  runs  (Fig.  13).  Turbulent 
luminous  streaks  are  evident  in  the  region  before  the  solenoid  where  the 
magnetic  field  strength  is  vanishingly  small.  As  the  pellet  enters  the  mag¬ 
netic  field  of  the  coil,  a  definite  change  in  the  structure  of  the  luminous  flow 
field  occurs.  The  luminous  streaks  which  characterize  the  turbulent  flow 
diminish  and  disappear  and  a  region  which  is  smooth  and  devoid  of  streaks 
appears.  The  pellet  position  at  which  the  maximum  magnetic  field  strength 
is  reached  is  shown  with  the  band  of  uncertainty.  The  lower  image  which 
runs  parallel  to  the  pellet  in  the  figure  was  usually  observed.  From  the 
information  obtained  in  this  first  set  of  runs,  it  appears  that  the  originally 
turbulent  wake  flow  is  substantially  affected  by  the  presence  of  the  mag¬ 
netic  field,  perhaps  suggestive  of  a  return  to  a  laminar  wake  flow. 

This  first  series  of  MHD  runs  raised  the  following  questions: 

A.  Was  the  disappearance  of  the  lumionous  streaks  as  the 
pellet  entered  the  coil  in  the  MHD  runs  possibly  related  1) 
to  an  end  effect  due  to  the  variation  of  the  magnetic  field  in¬ 
tensity  outside  of  a  solenoid;  2)  to  a  stagnation  point 
effect  in  which  the  heat  transfer  to  the  body  was  decreased 
by  MHD  effects,  thus  eliminating  the  luminous  "dye"  which 
previously  traced  the  flow  pattern;  3)  to  the  presence  of 
objects  in  the  wake  of  the  pellet-,  or  4)  to  some  other  effect, 
say  to  one  related  to  the  discharge  of  the  LC  circuit? 

B.  What  was  the  origin  of  the  second  luminous  streak  which 
runs  parallel  to  that  of  the  pellet? 

With  respect  to  question  A,  calculations  tend  to  indicate  that  the  mag¬ 
netic  field  end  effects  and  heat  transfer  effects  could  not  have  produced  the 
observed  large  changes  in  luminosity. 

A  second  set  of  experiments  was  made  to  help  answer  the  pre¬ 
ceding  questions.  The  diagnostic  instrumentation  was  improved  with  the 
addition  of  a  photoelectric  recorder  and  shadowgraph.  The  shadowgraph 
established  the  integrity  of  the  model  and  the  location  of  any  other  bodies 
relative  to  the  pellet.  Also,  a  technique  of  scribing  the  diaphragm  of  the 
light  gas  gun  was  developed  which  eliminated  all  debris  from  the  wake  of 
the  pellet.  With  respect  to  question  B,  additional  data  taken  with  streak 
camera  and  photoelectric  recorders  with  the  coil  unenergized  indicated 
that  the  second  image  is  real  and  correlates  with  the  focusing  of  the  re¬ 
flected  shock  wave  from  the  sides  of  the  pyrex  tube  test  section.  The  fo¬ 
cusing  reheats  the  wake  gas  producing  a  region  of  enhanced  luminosity 
which  trails  the  pellet  at  a  fixed  distance. 

The  effect  on  wake  luminosity  in  the  first  set  of  runs  shown  in  Figs.  12 
and  13  did  not  reappear  in  the  second  set.  A  sufficient  number  of  runs  were 
made  at  both  5.  0  and  2.  0  cm  Hg  pressure  to  verify  that  no  effect  was  evident¬ 
ly  present.  Essentially,  all  of  the  runs  in  the  second  set  with  and  without 
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magnetic  fields,  produced  drum  camera  photographs  which  appeared  as 
those  shown  in  Fig.  14. 

Two  significant  differences  existed  between  the  first  and  second  set 
of  runs.  First,  was  the  elimination  of  the  diaphragm  debris  from  the  wake 
of  the  pellet.  It  might  be  conjectured  that  an  interaction  of  the  debris  with 
the  magentic  field  in  some  way  caused  what  appeared  to  be  an  MHD  effect. 

The  diaphragm  is  0.  015  inches  thick  and  is  made  of  aluminum. 

Second,  there  was  an  uncontrollable  change  in  the  conditions  within 
the  range  between  the  first  and  second  set  of  runs.  The  first  set  of  runs 
were  made  while  a  sodium  vapor  oven  was  attached  just  upstream  of  the 
MHD  test  section.  This  was  due  to  the  fact  that  runs  for  a  sodium  schlieren 
experiments^  were  being  shot  alternately  with  the  base -MHD  runs.  Although 
the  oven  was  shut  off  for  the  base-MHD  experiment,  and  the  oven  cleaned, 
some  solid  sodium  may  still  have  remained  in  that  part  of  the  test  section. 
Between  the  first  and  second  set  of  the  base-MHD  experiment,  another 
ballistic  range  experiment  was  carried  out  which  required  a  clinically 
clean  test  section.  At  that  time,  the  sodium  vapor  oven  was  permanently 
removed  from  the  range,  and  the  range  was  thoroughly  cleaned. 

The  second  set  of  runs  for  the  base-MHD  experiment  took  place  in  a 
clean  range  in  contradistinction  to  the  first  set  which  took  place  in  a  range, 
a  part  of  which  may  have  been  coated  with  solid  sodium.  As  pointed  out  pre¬ 
viously,  the  assumption  of  easily  ionizable  particles  leading  to  Spitzer  con¬ 
ductivity  (the  largest  possible)  provides  Rjq  <  F.X,0>  while  the  lack 

of  such  particles  leads  to  conductivities  which  provides  Rjq  >  Rx,  0>  A-Z. 

Thus,  the  controlling  factor  may  have  been  the  respective  presence  and  then 
absence  of  sodium  dust  in  the  range. 

An  estimate  of  the  number  of  0.  1  micron  diameter  sodium  dust  par¬ 
ticles  required  to  produce  the  Spitzer  concentration  of  electrons  in  the  base 
region  has  been  made.  Particle  sizes  of  like  0.  1  micron  diameter  remain 
suspended  in  the  atmosphere  for  long  periods  of  time  due  to  the  effects  of 
Brownian  motion.  ^  A  sodium  ^here,  0.  1  micron  in  diameter  has  a  mass 
of  10~l6  grams  and  produces  10'  electrons  when  completely  ionized.  There¬ 
fore,  108  sodium  particles  per  cubic  centimeter  of  0.  1  micron  diameter  are 
required  to  produce  the  5  x  10^^  electrons  per  cubic  centimeter  for  Spitzer 
conductivity.  Engineering  handbook  values  for  average  dust  concentrations 
in  the  atmosphere  show  that  10^  particles  per  cubic  centimeter  is  well  with¬ 
in  engineering  expectations. 

Finally,  it  is  possible  that  (1)  one  of  the  other  assumptions  which 
went  into  the  calculations  is  unreasonably  optimistic  or,  (Z)  that  the  sug¬ 
gested  single  criterion  Rjq  <  Rx,0  is  too  simple  for  the  hypersonic  case. 
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APPENDIX  A 


MHD  Maps  for  Ballistic  Range  Experiment.  Argon  Only. 


The  MHD  parameters  in  the  base  region  of  a  pellet,  assuming  ionized 
argon  only,  are  calculated  using  the  following  assumptions; 

P,  =  5P 


h,  =  1/2  h  , 

b  '  stagnation 

Equilibrium  Chemistry. 


The  electrical  conductivity  of  the  ionized  argon  gas  is  calculated  using  the 
conductivity  equation  of  Lin,  Resler  and  Kantrowitz.  1  5  The  number  of  elec- 
J:rons  present  was  determined  from  the  degree  of  ionization  at  equilibrium 
from  Saha's  equation.  Th^’.  resultant  conductivity  from  Lin,  Resler  and 
Kantrowitz  differs  significantly  from  that  given  by  Spitzer  and  Harm  only 
when  the  temperature  is  below  6000°K. 

The  base  region  values  for  the  parameters  Hj^/Bjf,  and  Rjq  •  B 

are  shown  on  a  pressure  velocity  grid  in  Fig.  Al;  (S/^  )  ^  0»  =1.  0 

and  (RH)a)T=l.  0  sind  OJT /B  and  in  Fig.  A3. 
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Fig.  2  Schematic  of  the  transverse  MHD  vorticity  interaction  in  the  base 
region  of  a  sphere  with  flight  velocity  aligned  with  a  uniform 
magnetic  field. 
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Fig.  10  Laminar  Drum  Camera  Photograph.  Argon  Gas.  No  MHD. 
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Fig.  11  Turbulent  Drum  Camera  Photograph.  Argon  Gas.  No  MHD. 


-25- 


93  U 


FREE  STREAM  BASE  REGION  FREE  STREAM  BASE  REGION 


P«  *  5  cm  Hg 
V  *  13700  fps 

d  *  0.56  cm 


R  a  2700 

Rx 0-600 


Poo  =  5  c  m  H  g 
Px,o  =  lcm 

V  =  14260  fps 


B  =  1.26 


WEBER' 


M‘ 


R  «  2700 

Rh  *»  98 
Hm  »  30 


d  =  0.56  cm 


R|^  0.08 

WT  •«  0.8 


Fig.  12  Turbulent  Drum  Camera  Photograph.  Argon  Gas.  With  and 
Without  MHD.  FIRST  SET. 
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Laminar  and  Turbulent  Drum  Camera  Photographs  with  MHD. 
Argon  Gas.  SECOND  SET. 
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Fig.  A3  MHD  Base  Maps  -  Argon  Only. 
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